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Abstract. Neutrino oscillation is the only known phenomenon for physics beyond the standard 
model. To investigate this phenomenon, the understanding of low energy neutrino scattering (200 < 
E v < 2000 MeV) is the crucial task for high energy physicists. In this energy region, the charged 
current quasi-elastic (CCQE) neutrino interaction is the dominant process, and experiments require 
a precise model to predict signal samples. Using a high-statistics sample of muon neutrino CCQE 
events, MiniBooNE finds that a simple Fermi gas model, with appropriate adjustments, accurately 
characterizes the CCQE events on carbon. The extracted parameters include an effective axial mass, 
M A = 1.23 ± 0.20 GeV, and a Pauli-blocking parameter, K = 1.019 ± 0.01 1. 
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CCQE event selection in MiniBooNE 

The MiniBooNE detector, a spherical tank filled with mineral oil, detects Cerenkov 
light from charged particles Q. The identification of V^CCQE interactions relies solely 
on the detection of the primary muon Cerenkov light and the associated decay electron 
Cerenkov light in these events (Fig.Q]): 

v^+n^fi'+p ji ->e~ + V M + V e . 

By avoiding requirements on the outgoing proton kinematics, the selection is less de- 
pendent on nuclear models. The scintillation light from the proton, although not used 
directly in the v^CCQE analysis, is intensively studied in neutral current elastic scatter- 
ing events at MiniBooNE [3]. A total of 193,709 events pass the MiniBooNE v^CCQE 
selection criteria [2] from 5.58 x 10 protons on target collected between August 2002 
and December 2005. The cuts are estimated to be 35% efficient at selecting CCQE 
events in a 500 cm radius, with a CCQE purity of 74%. The predicted backgrounds are: 
74.8% CC ln+, 15.0% CC lrc , 4.0% NC \n ± , 2.6% CC multi-7z;, 0.9% NC elastic, 
0.8% v M CC \%~, 0.8% NC \k q , 0.6% r\lplK production, and 0.5% deep inelastic scat- 
tering and other events. Because pions can be absorbed via final state interactions in the 
target nucleus, a large fraction of the background events look like CCQE events in the 
MiniBooNE detector. "CCQE-like" events, all events with a muon and no pions in the 
final state, are predicted to be 84% of the sample after cuts. 



The detailed information of the Fermilab Booster neutrino beamline and the MiniBooNE neutrino 
detector are available elsewhere [Ql 0] . 



FIGURE 1. Schematic figure of a CCQE interaction. The primary Cerenkov light from the muon 
(Cerenkov 1 , first subevent) and subsequent Cerenkov light from the decayed electron (Cerenkov 2, second 
subevent) are used to tag the CCQE event. For most events, protons only emit scintillation light, and our 
selection is insensitive to this information. 

MiniBooNE CCQE events 

Fig [2] shows the data- Monte Carlo (MC) ratio of CCQE events as a function of muon 
kinetic energy T^(GeV) and muon scattering angle cosO^. Note the muon energy and 
muon scattering angle are the observables and the basis of all reconstructed kinematic 
variables in the v^CCQE channel. One can immediately see that the data-MC agreement 
is poor. There are 6 auxiliary lines: (a), (b), and (c) are equal neutrino energy lines, 0.4, 
0.8, and l.2GeV each, and (d), (e), and (f) are equal Q 2 lines, 0.2, 0.6, and l.OGeV 2 
each. The data-MC disagreement follows auxiliary lines of equal Q 2 , not equal neutrino 
energy, this indicates that the data-MC disagreement is not due to the neutrino flux 
prediction, but due to the neutrino interaction prediction, because the former is a function 
of neutrino energy and the latter is a function of Q 2 . So we assume that the data-MC 
disagreement comes from our neutrino interaction model and we adjust to the data. This 
is a critical task for MiniBooNE since the goal is to measure v e CCQE events, but MC 
and all reconstruction tools must be reliable and tested in copious v^CCQE events due 
to the blind analysis constraint on the v e CCQE channel. 

The data-MC disagreement is classified in 2 regions in this plane (Fig. [2]), 

1. data deficit at low Q 2 region, light gray band near the top left corner 

2. data excess at high Q 2 region, black band from the top right to the bottom left 

Since we are employing the Relativistic Fermi Gas (RFG) model [0] in our MC, we wish 
to fix these problems within the RFG model. The low Q 2 physics is usually controlled 
by nuclear model, so we want to tune the nuclear model, especially the strength of 
Pauli blocking, to fix region (1). This is justified because electron scattering data has 
not provided precise information for Pauli blocking in the low Q 2 region in terms of the 
RFG model For (2), we need to increase axial mass Ma to increase the cross section 
at high Q 2 . Here, the axial mass is understood as an experimental parameter in the axial 
form factor, 
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FIGURE 2. Ratio of MiniBooNE CCQE data/simulation as a function of reconstructed muon angle 
and kinetic energy. The prediction is prior to any CCQE model adjustments; the j 2 /dof = 79.5/53. The 
ratio forms a 2D surface whose values are represented by the gray scale, shown on the right. If the 
simulation modeled the data perfectly, the ratio would be unity everywhere. Contours of constant E v 
and Q 2 are overlaid. 



where gA is axial coupling constant (=1.267). This treatment is also justified because 
elastic electron scattering cannot measure the axial mass precisely. Interestingly, the high 
axial mass is also observed by the K2K experiment in Japan Q6Q . 



Pauli blocking parameter k 

Currently, MiniBooNE is using the NUANCE neutrino interaction generator |@]. In 
NUANCE, CCQE interactions on carbon are modeled by the Relativistic Fermi Gas 
(RFG) model Q4£| . To achieve our goal within the RFG model, we introduced a new 
parameter "kappa", K, 

E hl =^/p 2 +M 2 E 1o = k(^ p 2 +M 2 -co + E b ) (2) 

where M n is the target neutron mass, M p is the outgoing proton mass, Pp is Fermi 
momentum (=220MeV), Eg is binding energy (=34MeV), and (0 is the energy transfer. 
In the RFG model, is the energy of an initial nucleon on the Fermi surface and E\ 
is the lowest energy of an initial nucleon that leads to a final nucleon just above the 
Fermi momentum. The function of parameter K is to squeeze down the phase space of 
the nucleon Fermi sea, especially when the energy transfer is small. From Fig. [3j one 
can see that this parameter controls the Q 2 distribution only in the low Q 2 region. This 
is quite complementary to the role of Ma, since Ma mainly controls the Q 2 distribution 
in the high Q 2 region. 

We use these 2 parameters to perform a grid search to find the % 2 minimum. Here, we 
take into account all possible correlations between systematics by using the inverse of 
the full error matrix, not by adding systematics as pull terms. 
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FIGURE 3. Effect of and K variations on the MiniBooNE, flux integrated Q 2 distribution. The 
top plot shows various with fixed K, and bottom plot shows various K with fixed M^. Note, the 
variation has large impact at high Q 2 while the K variation has a significant impact only for Q 2 below 
~0.2GeV 2 . 



Fit result 

Finally, the parameters extracted from the MiniBooNE CCQE data are: 

M A = 1.23 ±0.20 GeV; (3) 
k = 1.019±0.011 . (4) 

Tab. CD shows the contributions to the systematic uncertainties on Ma and K. The 
detector model uncertainties dominate the error in Ma due to their impact on the energy 
and angular reconstruction of CCQE events in the MiniBooNE detector. The dominant 
error on K is the uncertainty in the Q 2 shape of background events. 

The result of this fitting, including all sources of systematic uncertainty, is shown in 
Fig. SI Since the background error dominates at low Q 2 , and it drives the large error 
bars at low Q 2 . Note that, the shape uncertainty of the background, namely the Q 2 
distribution shape uncertainty of CCln + events, is not included in these error bands. 
From the data, we know that the predicted Q 2 shape of CC l7T + events have large 
errors The extracted shape information from the data is implemented in our MC, 
and the fit is performed again. The result of those 2 fits, one using the MC predicted 
CCln + distribution, and the other is using MC tuned on CCln + data, are shown with 
the star and the triangle in the inserted plot in Fig. HI The difference is interpreted as a 
background shape uncertainty error and added to the extracted parameters. 

Fig. [5] shows the agreement between data and simulation after incorporating the Ma 
and K values from the Q 2 fit to MiniBooNE v^CCQE data. Comparing to Fig. [21 the 
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FIGURE 4. Reconstructed Q 2 for v^CCQE events including systematic errors. The simulation, before 
(dashed) and after (solid) the fit, is normalized to data. The dotted (dot-dash) curve shows backgrounds 
that are not CCQE (not "CCQE-like"). The inset shows the Id CL contour for the best-fit parameters 
(star), along with the starting values (circle), and fit results after varying the background shape (triangle). 

TABLE 1. Uncertainties in and K from 
the fit to MiniBooNE v M CCQE data. The total 
error is not a simple quadrature sum because 
of the correlation between the two parameters. 



error source 


8M A 


5k 


data statistics 


0.03 


0.003 


neutrino flux 


0.04 


0.003 


neutrino cross sections 


0.06 


0.004 


detector model 


0.10 


0.003 


CC 7Z + background shape 


0.02 


0.007 


total error 


0.20 


0.011 



improvement is substantial and the data are well-described throughout the kinematic 
phase space. Since the whole kinematic space is fixed, not surprisingly, all of the 
individual kinematic variables exhibit good data-MC agreement. Fig. [6] shows that data 
and MC agree well within error bars for reconstructed muon neutrino energy and muon 
scattering angle. 

In general, varying Ma allows us to reproduce the high Q 2 behavior of the ob- 
served data events. A fit for Ma above Q 2 > 0.25 GeV 2 yields consistent results, 
Ma = 1.25 ±0.12 GeV (Fig. [7]). However, fits varying only Ma across the entire Q 2 
range leave considerable disagreement at low Q 2 . This data-MC disagreement at low 
Q 2 would eventually reflect in data-MC disagreement in reconstructed neutrino energy, 
because data-MC disagreement in Q 2 spreads out in the kinematic plane and would af- 
fect the energy reconstruction across a wide region. The Pauli-blocking parameter K is 
instrumental here, enabling this model to match the behavior of the data down to Q 2 = 0. 
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FIGURE 5. Ratio of data/simulation as a function of muon kinetic energy and angle after the CCQE 
model adj ustments ; the ^ 2 /dof = 45.1/53. Compare to Figure [2] 




FIGURE 6. Reconstructed muon neutrino energy and measured muon scattering angle, line notations 
are the same as Fig [4] 
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FIGURE 7. Reconstructed Q 2 for v M CCQE events and data-MC ratio in the kinematic plane. The left 
plot is the analogy of Fig. [4] and the right plot is the analogy of Fig. |5]but the fit is performed using 
only, with fixed K (=1.0, no enhanced Pauli blocking). 
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FIGURE 8. Reconstructed Q 2 for v^CCQE events, including statistics error only. The simulation, old 
model (dashed) and new model (solid) is normalized to data. 



Anti-neutrino CCQE preliminary result 

Finally, we tested the modified RFG model in a new sample of MiniBooNE antineu- 
trino data (for other preliminary results for antineutrino run, see [|9D). If our assumption 
is correct, this RFG model should also succeed in modeling v^CCQE events. The result 
is shown in Fig. [8l Although statistics are low, one can tell the new model will describe 
the features of this data better than the original model. 

In summary, taking advantage of the high-statistics MiniBooNE CCQE data, we 
have extracted values of an effective axial mass parameter, Ma = 1.23 ±0.20 GeV, 
and a Pauli-blocking parameter, K = 1.019 ±0.011, achieving substantially improved 
agreement with the observed kinematic distributions in this data set. 

The Ma value reported here should be considered an "effective parameter" in the sense 
that it may be incorporating nuclear effects not otherwise included in the RFG model. 
Future efforts will explore how the value of Ma extracted from the MiniBooNE data is 



altered upon replacement of the RFG model with more advanced nuclear models [10]. 
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